Charm and bottom production in inclusive double Pomeron exchange in heavy ion 

collisions at the LHC 



O 

O 

Q 



43 ■ 

Oh 1 
i ' 

43 . 



> 

o 

oo 

o 
o 



X 



M. B. Gay Ducati, M. M. Machado and M. V. T. Machado 
High Energy Physics Phenomenology Group, GFPAE, IF-UFRGS 
Caixa Postal 15051, CEP 91501-970, Porto Alegre, RS, Brazil 

The inclusive double Pomeron exchange cross section for heavy quark pair production is calcu- 
lated for nucleus-nucleus collisions at the LHC. The present estimate is based on hard diffractive 
factorization, corrected by absorptive corrections and nuclear effects. The theoretical uncertain- 
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production channels giving a similar final state configuration are discussed as well. 
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I. INTRODUCTION 

At the LHC, heavy quark pairs are produced in large 
quantities and they are important both for perturba- 
tive QCD physics study and for understanding the back- 
ground to other processes. For instance, heavy flavored 
hadrons may also produce high momentum leptons, and 
therefore they always constitute a potential background 
to new physics. The theory of hadronic production of 
heavy quarks is in reasonably good shape. The inclu- 
sive production is known beyond the next-to-leading level 
order in QCD and the results are considered well estab- 
lished. Recently, the clean topologies of exclusive particle 
production in electromagnetic interactions in context of 
hadron-hadron and nucleus-nucleus collisions mediated 
by colorless exchanges such the QCD Pomeron have at- 
tracted an increasing interest [1]. The cross sections for 
these processes are smaller than the correspondent inclu- 
sive production channels, which it is compensated by a 
more favorable signal/background relation. Experimen- 
tally, exclusive events are identified by large rapidity gaps 
on both sides of the produced central system and the 
survival of both initial state particles scattered at very 
forward angles with respect to the beam. 

Our goal in this work is to estimate the inclusive double 
Pomeron exchange cross section for heavy quarks pro- 
duction in heavy ion collisions at the LHC. We focus 
on the coherent diffractive production of heavy quarks 
in relativistic nucleus- nucleus collisions, i.e. the process 
A + A -> A + [LRG] + [QQ + X] + [LEG] + A, where 
[LRG] stands for Large Rapidity Gap. Notice that the 
current study is complementary to our previous investi- 
gations on single diffractive production of charm/bottom 
presented in Ref. 0. It is also possible to occur inco- 
herent diffractive scattering, where one or both nuclei 
will be excited, that subsequently decays into a system 
of colorless protons, neutrons and nuclei debris. Such 
processes are of experimental interest at the LHC, for 
instance at ALICE experiment [3|]. Here, we start by 
the hard diffractive factorization, where the diffractive 
cross section is the convolution of the diffractive parton 
distribution functions and the corresponding diffractive 
coefficient functions, in a similar way as for the inclusive 



case. To further correct this approach it is introduced the 
appropriate absorptive (unitarity) effects which cause the 
suppression of any large rapidity gap process being im- 
portant for the reliability of predictions It is also 
considered the nuclear effects for nucleus-nucleus colli- 
sions and the predictions compared to other approaches 
available in the literature. 

This paper is organized as follows. In the next section 
it is summarized the main formulas considered to com- 
pute the inclusive double Pomeron exchange (DPE) cross 
section and diffractive ratios for the hadroproduction of 
charm and bottom. We also present the procedure to 
evaluate the nuclear coherent diffractive cross sections at 
the LHC energies. At this point it is introduced other 
models for the inclusive diffractive heavy quark produc- 
tion. In the last section, the numerical results are shown 
for the inclusive and diffractive cross sections as a func- 
tion of energy and give predictions for the corresponding 
diffractive ratios for pp and AA collisions. Discussion 
on the nuclear dependence of cross sections and corre- 
sponding suppression factors are addressed. A compari- 
son to competing channel, as the exclusive photo-nuclear 
reactions in AA collisions and two-photon production, is 
performed. 



II. HEAVY-QUARK PRODUCTION IN 
INCLUSIVE DPE PROCESS 



Let us present the main formulas for the inclusive 
diffractive cross sections for the production of heavy 
quarks in proton-proton collisions at high energies. 
We rely on diffractive factorization formalism Q and 
take into account both the absorption effects (multiple 
Pomeron exchange) and nuclear effects using Glauber 
based approach. The limitations and theoretical in- 
completeness of such a picture are well known. How- 
ever, it is reasonable for a first exploratory study. In 
the inclusive case, the process is described for partons 
of two protons, interacting to produce a heavy quark 
pair, p + p — > QQ + X , with center of mass energy -y/5. 
At LHC energies, the gluon fusion channel dominates 
over the qq annihilation process and qg scattering. The 
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NLO cross section is obtained by convoluting the par- 
tonic cross section with the parton distribution function 
(PDF), g(x, jtxj?), in the proton, where pf is the factoriza- 
tion scale. At any order, the partonic cross section may 
be expressed in terms of dimensionless scaling functions 
fh l that depend only on the variable p 0, 
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hJ = Q,q,9, specifying the types of 
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the annihilating partons, s is the partonic center of mass, 
rriQ is the heavy quark mass, pr is the renormalization 
scale (sq — 1 GeV 2 ). It is calculated as an expansion 
in powers of a s with k = corresponding to the Born 
cross section at order 0{a 2 s ). The first correction, k — 1, 
corresponds to the NLO cross section at O(al). To calcu- 
late the fij in perturbation theory, both renormalisation 
and factorisation scale of mass singularities must be per- 
formed. The subtractions required are done at the mass 
scale p. The running of the coupling constant a s is deter- 
mined by the renormalization group. The total hadronic 
cross section for the heavy quark production is obtained 
by convoluting the total partonic cross section with the 
parton distribution functions of the initial hadrons Q 



(Tpp(s,TOQ) 
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dx 2 f?(x 1 ,p F )fP(x 2 ,p F ) 



x a 2J (s,m Q ,p F ,p R ), 



(2) 



with the sum i, j over all massless partons. Here, x\ y2 
are the hadron momentum fractions carried by the inter- 
acting partons, f?^ is the parton distribution functions, 
evaluated at the factorization scale and assumed to be 
equal to the renormalization scale in our calculations. 
Here, the cross sections were calculated with the follow- 
ing mass and scale parameters: p c = 2m c , m c = 1.5 
GeV, pb = mb = 4.5 GeV, based on the current phe- 
nomenology for heavy quark hadroproduction Q . 

For diffractive processes, we rely on the hard diffractive 
factorization [Ej], where the Pomeron structure (quark 
and gluon content) is probed. In the case of single 
diffraction, a Pomeron is emitted by one of the colliding 
hadrons. That hadron is detected, at least in principle, in 
the final state and the remaining hadron scatters off the 
emitted Pomeron. A typical single diffractive reaction is 
given by p + p — > p + QQ + X, with the cross section as- 
sumed to factorise into the total Pomeron-hadron cross 
section and the Pomeron flux factor [5], fw/i(xw\ I^D- 
As usual, the Pomeron kinematical variable is defined 



as ip = s^ 1 /sij, where 



,0) 



is the center-of-mass en- 



ergy in the Pomeron-hadron j system and ^fsu = \fs the 



center-of-mass energy in the hadron i-hadron j system. 



The momentum transfer in the hadron i vertex is denoted 
by ij. A similar approach can also be applied to double 
Pomeron exchange (DPE) process, where both colliding 
hadrons can in principle be detected in the final state. 
Thus, a typical reaction would be p+p — > p+QQ+X +p, 
and DPE events are characterized by two quasi-elastic 
hadrons with rapidity gaps between them and the central 
heavy flavor products. The inclusive DPE cross section 
may then be written as, 



da(pp -> pp + QQ + X) 
dx$ dx { $ d\h\d\t 2 \ 



fTP/p{x^\ |*l|) fw/p&w', |*2 j) 



(2) 



x a (TP + TP ^QQ + X), 

i,3=<l,9 

(3) 

where the Pomeron-Pomeron cross section is given by, 

a (TP + TP ^ QQ + X) = J J dx 1 dx 2 a. lJ {s,m 2 Q ,p 2 ) 

x Si/sp (Pi,p 2 ) f j/1P (/? 2 ,M 2 ) ,(4) 

where fi/Tp(/3,p 2 ) are the diffractive parton (quark, 
gluon) distribution functions (DPDFs) evaluated for par- 
ton momentum fraction f5 a — x a /x% (a — 1,2) and evo- 
lution scale p 2 . 

We further correct Eq. ([3]) by considering the sup- 
pression of the hard diffractive cross section by multiplc- 
Pomeron scattering effects (absorptive corrections). This 
is taken into account through a gap survival probabil- 
ity, <Sg ap , which can be described in terms of screening 
or absorptive corrections [9(. There are intense theo- 
retical investigations on this subject in the last years 
and we quote Ref. [l(| for a discussion on several 
theoretical estimations for the gap survival probabili- 
ties. As a baseline value, we follow Ref. [11] . For 
the present purpose, we consider S 2 ap = 0.032 (0.031) at 
y/s = 5.5 (6.3) TeV in nucleon-nucleon collisions, which 
is obtained using a parametric interpolation formula for 
the KMR survival probability factor [ll| in the form 
•Sgap = a/[b+\n(-y/s/s )} with a = 0.126, b = -4.688 and 
,sq = 1 GeV 2 . This formula interpolates between survival 
probabilities for central diffraction (CD) in proton-proton 
collisions of 4.5 % at Tevatron and 2.6 % at the LHC. 

For the heavy quarks production in nucleus-nucleus 
collisions, two processes can occur: (a) coherent diffrac- 
tive scattering, A+A —> A+QQ+X +A, and (b) incoher- 
ent diffractive scattering, A + A ->• A* + QQ + X + A(*) . 
In the first case, both nuclei emit Pomerons and par- 
tons from them interact with each other. Thus, both 
nuclei remain intact at the final state. For the incoher- 
ent diffractive scattering, one or both nuclei can be ex- 
cited, represented by A*, that subsequently decay into a 
system of colorless protons, neutrons and nuclei debris. 
Here, we focus on the coherent case. In order to calculate 
the AA cross section the procedure presented originally 
in Ref. [l^] (see Refs. [l3| for further applications to 
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soft processes), the so-called criterion C, will be used. 
The central diffraction (coherent) cross section for AB 
collisions is given by: 



CD 
'AB 



a AB \ a PV ) 



271 t % 

a AB \ p 



,CD\ 



(5) 



where a l AB is the inelastic AB cross section considered 
as a function of the nucleon-nucleon total cross section a 
(cr™ and (jpp are the inelastic and CD cross sections 
in proton-proton case, respectively). The application 
of the formalism to AA collisions turns out to be dif- 
ficult due to the absence of an explicit expression for 
°'aa(' 7 )- I n R- er - [HL1IH the authors considered optical 
approximation in which a l AA is given by corresponding 
formula for pA with A — > AB and an effective profile 
function for two colliding nuclei at impact parameter b, 
T AB = J d 2 bT A (b)T B (b-b). The final expression for CD 
cross section in AA collisions is given by 



°aa E = A 2 J d 2 bT AA (b) exp [-A 2 a™ T AA (b)] a 



DPE 

PP ■ 



(0) 



where we consider Woods-Saxon nuclear densities and 
the inelastic cross section at the LHC is taken to be 
a™ (VS = 6 TeV) = 73 mb [H El. 

To analyze the model dependence of the cross sec- 
tion, we consider another approach to inclusive diffrac- 
tive production of heavy quarks. In order to do so, the 
Bialas-Landshoff (BL) approach [lj], HH for the process 
p + p — > p + QQ + p is taken into account. The calcula- 
tion that follows concerns central inclusive process, where 
the QCD radiation accompanying the produced object is 
allowed. Thus, we did not include a Sudakov survival 
factor T(k,/i) [16[ which is needed for exclusive central 
processes. The cross-section is given by (l7| : 

(T]p]p(BL) = — — !— s fWJ^[F(h)F(t 2 )] 2 dPH, (7) 
Is (2ir) J 

where F (t) is the nucleon form-factor approximated by 
F (t) — exp(bt), with slope parameter b — 2 GeV~ 2 . 
The differential phase-space factor dPH has the form, 

dPH = d 4 kx5 (kf) d 4 k 2 S (fcf ) d 4 n S (rf - m 2 Q ) 

x d 4 r 2 5 (r 2 - m%) 9 (A?) 6 (*°) 6 (r?) 6 (r°) 



<5 (4) (pi + p 2 - h - k 2 - n - r 2 ) 



(8) 



where tuq is the mass of the produced quarks. Following 
[l7| . the use of Sudakov parameterization for momenta 
is given by 



x y 
Q = -Pi + -P2 + v, 

s s 

k 2 = — p\ + y 2 p 2 + v 2 , 
s 



ki = xipi H p 2 + Vi, 

s 

T2 = XQPl +VQP 2 +VQ, 



where v, V\, v 2l vq are two-dimensional four- vectors de- 
scribing the transverse components of the momenta. The 
momenta for the incoming (outgoing) protons are p\ , p 2 



QQ 


(Tine [fJ-b] 


odpe \pb] -Rdpe [%] 


cc 


7811 


13.6-0.53 0.17-7xl0 -3 


bb 


393 


0.053-0.027 0.01-0.007 



TABLE I: The inclusive and DPE (corrected by absorption 
effects) cross sections in pp collisions at the LHC. For the 
inclusive diffractive cross section the first value corresponds 
to the Ingelman-Schlein approach and the second one the 
Bialas-Landshoff approach. The corresponding diffractive ra- 
tios, -Rdpe, are also presented. 



(fci, k 2 ) and the momentum for the produced quark (an- 
tiquark) is r 2 [r\ ) , whereas the momentum for one of the 
exchanged gluons is Q. The square of the invariant ma- 
trix element averaged over initial spins and summed over 
final spins is given by (l7| . 



\M fi 



xiy 2 H 



{sxqVq) \0id 2 ) b 1 x b 2 
( 4 m-k \ 

exp [2(3 (h + t 2 )] 



s8\8 2 



(9) 



In the expression above, <5i — 1 — xi, 5 2 — 1 — y 2 , t\ = 
—v 2 and t 2 = —v\. The factor exp [2(3 (t± + t 2 )] takes into 
account the effect of the momentum transfer dependence 
of the non-perturbative gluon propagator with (3 = 1 
GeV -2 . The overall normalization can be expressed as, 



H = S 2 ap x 2s 



4Trm Q (G 2 D ) 3 M 4 



-, 2 



9 (2tt) 2 



(10) 



where a s is the perturbative coupling constant (it de- 
pends on the hard scale) and ag (supposed to be in- 
dependent of the hard scale) is the unknown nonper- 
turbative coupling constant. In the numerical calcula- 
tion, we use the parameters [l?} e = 0.08, a' = 0.25 
GeV- 2 , n = 1.1 GeV and G 2 D = 30 GeV-V 1 - The 
Regge Pomeron trajectory is then ajp(i) = 1 + e + a't. 
It is taken fc min = for the minimum value for the 
transverse momentum of the quark. For the strong cou- 
pling constant, we use a s = 0.2 (0.17) for charm (bot- 
tom). An indirect determination of the unknown param- 
eter ao has been found in Ref. [l8j using experimental 
data for central inclusive dijet production cross section 
at Tevatron. Namely, it has been found the constraint 
S 2 ap (-y/s = 2TeV)/a 2 = 0.6, where S 2 ap is the gap sur- 
vival probability factor (absorption factor). Considering 
the KMR [TO] value 5^ ap = 0.045 for CD processes at 
Tevatron energy, one obtains a 2 , = 0.075. For the nu- 
clear version of the cross section we use the same proce- 
dure described before. 

In what follows, we present the numerical results using 
procedures referred above and compare them to distinct 
theoretical approaches in literature. Competing channels 
producing the same final state configuration (two large 
rapidity gaps and central produced system) are also dis- 
cussed. 
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A. Results and discussion 

First, we present the results for the inclusive and 
diffractive heavy quarks cross sections for hadronic col- 
lisions at LHC. The calculations for the inclusive and 
diffractive cross sections as well as the diffractive ratios 
to heavy quark production in proton-proton collisions are 
showed at Tab. (|l]). For the inclusive diffractive cross 
section the first value corresponds to the partonic pic- 
ture of Pomeron, Eq. <j3j) , and the second one to the 
Bialas-Landshoff approach, Eq. (0. We assume the 
value Sg ap — 0.026 for the absorption corrections at en- 
ergy of 14 TeV. The partonic PDFs and scales are men- 
tioned in previous section. For the diffractive gluon PDF, 
we take the experimental (HI collaboration) FIT A (l9j . 
The main theoretical uncertainty in the diffractive ratio 
is the survival probability factor, whereas uncertainties 
associated to factorization/renormalization scale, parton 
PDFs and quark mass are minimized taking a ratio. The 
present results are consistent with a previous estimate 
performed in Ref. [20j . where a value 5g ap = 0.04 was 
considered and cross sections were computed at LO ac- 
curacy. 

It is important to discuss the uncertainties in our es- 
timates presented in Tables I and II. In the partonic 
Pomeron model, the main theoretical uncertainties come 
from the factorization/renormalization scale and from 
the diffractive PDFs. The scale dependence is stronger 
for charm case (the error band reaches a factor around 
2 (H) and very stable for bottom. We have checked that 
the cross sections are insensitive to a distinct choice on 
the diffractive gluon PDF (for instance, HI Collabora- 
tion, FIT B). The main theoretical uncertainty in the 
Bialas-Landshoff model is the nonperturbative coupling 
ato- We have used the adjusted value from [l8|, extracted 
from experimental data for central inclusive dijet cross 
section at y/s = 1.8 TeV and we believe that it is sta- 
ble for extrapolations to high energies. A considerable 
concerning in our calculation is the assumption of the 
gap survival probability in AA collisions being the the 
same as for pp collisions and in addition to be final- 
state independent. The last point is better understood 
than the first one. Considering for simplicity the usual 
eikonal (one channel) model [21|, it is known that the 
gap survival probability for DPE processes is given by 
S 2 ap = a[2i/(s)]~ a j(a, 2v(s)), where ~f(a,2v) denotes the 
incomplete Euler gamma function. The variable a de- 
pends on energy, s, and on final state mass, M. Namely, 
a(s,M) = 2R 2 {s)/R 2 {^ I ) 1 with R = 4[JJ§ +a' ln(s/s )]. 
In addition, v(s) ~ 2 -ni{ i (s) ( s / s oY ■ It is clear that 
a(s, M) has a rather weak dependence on M 2 as it is 
proportional to a' ln(A/ 2 ) (with a' small) over a rela- 
tively narrow domain. Therefore, it can be factored out 
of M 2 integration. For a numerical example of the small 
sensitivity of final-state configuration, in [2l| the gap sur- 
vival probability is computed for hard diffractive dijets 
and for meson x c (3415) and the deviation is shown to 
be of order 5% at 14 TeV. The question about the gap 



survival in heavy-ion collisions is an open question in lit- 
erature. It is naively expected that it depends on impact 
parameter and will be different in pp and AA collisions. 
For the pp case, it has been shown in Ref. [22j that the 
uncertainties on the numerical predictions for S 2 ap come 
from the transverse spatial distribution of gluons enter- 
ing in Phard(i>). This quantity is the probability for two 
gluons to collide at same transverse point as a function of 
the pp impact parameter, given by the convolution of the 
transverse spatial distributions of the gluons in the col- 
liding protons. An additional source of uncertainty arises 
from the modeling of the pp elastic amplitude, T(b). It 
was also shown in [22[ that moderate extra suppression 
results from fluctuations of the partonic configurations of 
the colliding protons and at LHC energies absorptive in- 
teractions of hard spectator partons associated with the 
process of gluon fusion producing a massive final state 
reach to the black-disc regime and it would generates 
sizable additional suppression. A careful analysis along 
these line was not still done for A A collisions. 

Let us discuss the model dependence of the cross sec- 
tions. The partonic approach for Pomeron produces a 
larger cross section compared to Bialas-Landshoff model. 
This comes mostly from the distinct energy dependence 
in the phenomenological models. The first one considers 
a semihard Pomeron intercept fitted from HERA diffrac- 
tive DIS data, where e ~ 0.12. On the other hand, 
the soft Pomeron models consider the conservative soft 
Pomeron intercept, e = 0.08. The deviation less evident 
for bottom may come from a distinct toq dependence on 
each approach. The present calculation can be compared 
to the exclusive diffractive (central diffraction, CD) chan- 
nel, p + p — > p + QQ + p. Very recently, in Refs. [23| it 
has been computed considering the KKMR procedure to 
compute the CD processes for open charm and bottom 
production. Notice that the Ingelman-Schlein approach 
(IS) gives larger cross sections so far as it refers to the 
inclusive central diffraction, p+p —> p + QQX +p, where 
the remnants in the final state of the Pomeron-Pomeron 
sub-system (X) would share the energy and momentum 
continuously, and therefore no distinct gap would emerge. 

As a cross-check for Tevatron at y/s = 1.96 TeV we got 
1.08-0.37 /ib for charm and 6-16 nb for bottom. These 
values can be compared to the exclusive heavy quark 
production recently computed using the KKMR proce- 
dure 23[. In that case, the imaginary part of the ampli- 
tude of the exclusive diffractive quark pair production, 
p + p -» p + QQ + p, is given by 23] , 



M 



s n 2 S c 



QA « 2{N 2 - 1) 



d kq- 



V 



A Q A Q 



„2 2 2 



X J 7 g (x 1 ,x' 1 ,K ,K 1 ,t 1 ) J r g(x 2 ,x' 2 ,K ,K 2 ,t 2 ), (11) 

where Xq and Xq are helicities of heavy Q and Q, re- 
spectively. Above, the quantities T g are the off-diagonal 
unintegrated gluon distributions (UGDs) in each nucleon, 
K o,i,2 x 2) are the gluon transverse momenta (longitu- 
dinal momentum fractions of active gluons) in the KKMR 
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CaCa [cc] 


PbPb [cc] 


CaCa [66] 


PbPb [66] 




22.8-2.8 


31.1-4.2 


0.25-0.14 


0.32-0.2 




3 - 0.4 x 1CT 4 


2 - 0.2 x 10~ 4 


8 - 4 x 10~ 5 


4-3 x 10~ 6 



TABLE II: The diffractive coherent cross section for CaCa and PbPb collisions. The first value corresponds to the partonic 
Pomeron approach and the second one to the Bialas-Landshoff approach. Diffractive ratios are also presented. 



approach and i\ t2 are the quark (antiquark) momenta. 
The vertex factor V£^\_(ki, K2,ii,£ 2 ) is the production 
amplitude of a pair of massive quarks with helicities Xq q 
and momenta There is a lot of uncertainties con- 

cerning the off-diagonal UGDs even in the nucleon case. 
After some exercise, we are able to estimate that in that 
work the cross section is ss 2.97-1.1 fib for charm case 
(at Tevatron) and sa 1.3 nb for bottom at the LHC. 

Let us move now to the computation of nucleus-nucleus 
cross sections. The results are shown in Tab. Ql] for 
diffractive coherent collisions. The corresponding diffrac- 
tive ratio is obtained from 



T DPE 

AA 



_mc 
AA 



100%, 



(12) 



where (Jaa is given by the hadroproduction cross section 
in Eq. (|2|) and using minimum bias cross section for AA 
collisions. The diffractive DPE cross section in nucleus- 
nucleus collisions is obtained from Eq. 

In Tab. ((ill, the coherent cross section is presented 
for calcium and lead nuclei. The general trend on their 
magnitudes follows the main features occurring at the 
proton-proton case. Namely, larger deviation between IS 
and BL models for the charm case, whereas the difference 
diminishes for bottom. The diffractive ratios are quite 
small compared to the proton case. The present calcula- 
tion is still somewhat consistent with earlier calculation 
presented in Ref. 24] (absorption is not take into ac- 
count in that work), which predicts cr DPE = 0.02 — 60 
fib for charm and <r DPE = 0.003 — 1.4 fib for bottom 
in PbPb collisions. The nuclear effect is proportional 
to A 1 / 3 , which enhances by one order of magnitude the 
cross section compared to the nucleon case. The inco- 
herent case is not addressed here, but a procedure to 
its estimation has been presented in our previous study 
on single diffractive quark pair production Q. As the 
nuclear shadowing correction is an important theoretical 
uncertainty in the current work, it is timely to compare 
the predictions with different models. One possible al- 
ternate approach is the formalism considered by Muller 
and Schramm 25j for the meson (and Higgs) production 
in double Pomeron exchange. In that case, the cross sec- 
tion for particle production via two Pomerons exchange 
is written as 



a AA = J f dxi dx 2 fv(xi) f w (x 2 ) (Twtp(stpw), 



(13) 



where fw(x) is the distribution function that describes 
the probability of finding a Pomeron in the nucleus with 



energy fraction x and crpp^jpip) is the subprocess cross 
section with squared energy, Spjp. It has been shown in 
[25| that function /p, integrated over momentum transfer 
t, can be written as 



3Af3 Q 
2tt 



exp 



x 2 M\ 

Ql 



(14) 



with the soft Pomeron parameters e = 0.085 and (3q = 1.8 
GeV -1 . Assuming a Gaussian expression for the nuclear 
form factor in impact parameter space, one has Qo = 60 
MeV. In addition, m p and Ma are the proton and nu- 
cleus masses, respectively. The differential cross section 
in impact parameter space is then given by, 



d^r A 



d 2 b 



\ 2n z J J J xi x 2 



x exp 



M A ( 2 

Ql [Xl 



4) 



exp 



b 2 Ql 



X\X 2 S' 



2< 



]PIP— >X ( \ 

a AA \x x x 2 s) 



(15) 



To compute the cross section for the subprocess IP IP — > X, 
authors of Ref. [25[ rely on the Pomeron model of 
Donnachie-Landshoff. Namely, it is assumed that the 
Pomeron couples to the quarks like an isoscalar photon. 
This means that the cross section of such a subprocess 
can be obtained from suitable modifications on the cross 
section for 77 — > X. Another aspect to be considered 
is that the Pomeron-quark-quark vertex is not point-like, 
and when either or both of the two quark legs in this ver- 
tex goes far off shell, the coupling is known to decrease. 
So the quark-Pomeron coupling /?o must be replaced by 

MQ 2 ) = A^o/^o + Q 2 )> wher e Mo = 1-2 GeV 2 is a 
mass scale characteristic of the Pomeron; in our case of 
heavy quark production we have set Q = tuq . Therefore, 
the process IP IP — > X is totally similar to the one initi- 
ated by photons unless from an appropriate change of 
factors. The cross section c(lPIP — > QQ) is obtained by 
considering the cross section (7(77 — > QQ) by changing 

the fine-structure constant squared, a 2 , by 9 /3q / (167r 2 ). 
We consider the direct QED contribution to the process 
77 ~* QQ-, which should be dominant contribution in the 
present purpose. The corresponding DPE cross section 
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Heavy Quark JPW (IS) JPW (BL) 



77 



cc 
bb 



32.5 fib 
0.32 /ib 



4.2 (ib 
0.2 /xb 



1.88 mb 59 mb 
2.1 /ib 10 (ib 



TABLE III: Comparison of total cross sections for DPE, two- 
photon and photon-Pomeron channels for central charm and 
bottom production in PbPb collisions at the LHC. 



is given by 
a(WTP -> QQ) 



= 00 



Ami 



In- 



1 -£J 



U 1 + 



(16) 



where cr = 27eg/3 4 /(47r) and cj = wl — (4mg/s). In- 
tegrating over impact parameter we found the following 
cross sections: 0.12 mb for charm and 0.8 nb for bot- 
tom in PbPb collisions at the LHC. This is consistent 
with the computed cross sections for DPE production 
of 77-mesons in Ref. which gives 0.76 mb and 0.81 

nb for rj c and rjb, respectively. After correcting them 
by a gap survival probability one obtains ct^(cc) = 3.8 
(ib and a^(bb) = 0.02 nb. The strong suppression for 
bottom case comes from the model for the Pomeron- 
quark coupling, producing a suppression proportional to 
(/ip + m^) 4 /(/iQ + ml) 4 ~ 7 x 10~ 4 plus the suppres- 
sion due to distinct quark charge e 4 /e 4 = 1/16. For the 
charm case, the Muller and Schramm model produces a 
similar result as for the Bialas-Landshoff model, whereas 
produces a smaller cross section for bottom case. 

Finally, we compare the present calculation with simi- 
lar channels of production. Let us consider first the exclu- 
sive photoproduction of heavy quarks in ultra-peripheral 
collisions in heavy ion collisions [l[, which is dominated 
by photon-Pomeron interactions. The reason is that the 
final state configuration is similar to the coherent diffrac- 
tion (both nuclei remain intact and one rapidity gap). 
It has been recently computed in Ref. [27j for proton- 
proton collisions. It was found a^ c = 0.16 — 0.53 /ib 
for charm and a^ c — 0.32 — 3 nb for bottom, where the 
band correspond to model dependence using distinct sat- 
uration models. These values can be compared to Table 
I and it is verified that the diffractive channel dominates 
over exclusive photoproduction channel in hadronic col- 
lisions. On the other hand, for heavy ion collisions it 
was found in Ref. [28| that Cp£p b = 59 mb for charm 
and crp^p b = 10 (ib for bottom (lower bound cross sec- 
tion relying on saturation physics). Therefore, for heavy 
ions the photon-Pomeron channel should dominate over 
the diffractive channel (see Table II). Here, some com- 
ments are in order. In general, the parton saturation 
models underestimate the experimental results for heavy 
quark photoproduction at DESY-HERA regime. The 
reason is that they present a smaller effective power on 
energy when compared to models including QCD evolu- 
tion in gluon distribution. Therefore, the results com- 



ing from saturation physics can be considered a lower 
bound for the ultra-peripheral cross sections. The en- 
hancement of UPC's compared to single diffractive chan- 
nel in AA reactions is easily understood from their dis- 
tinct A-dependences: photoproduction grows as oc .4 4 
(Z 2 enhancement from the equivalent photon flux plus 
an additional enhancement from the photonuclear cross 
section, c^f °c A 4 / 3 ), whereas single diffractive cross 
section gets a factor A 2 S ~ A 1 / 3 . Another channel pro- 
ducing similar final state configuration is the two-photon 
channel. For instance, it has been calculated in Ref. [29[ 
using the color dipole formalism for the photon-photon 
interaction. The numerical results for these calculations 
are cr^ pb = 1.81-1.95 mb for charm and a^ ph = 2-2.2 
(ib for bottom, using m c = 1.3 GeV and mb — 4.5 GeV. 
The lower value correspond to a parton saturation model 
and the upper one to the BFKL approach. The average 
value for the 2-photon cross section is presented in Table 
III. Thus, for heavy ions collisions also the two-photon 
induced processes are larger than the diffractive channel. 
For experimental separation of distinct channels, further 
kinematical cuts have to be imposed. As a summary of 
distinct production channels, in Table III the cross sec- 
tions are presented for PbPb collisions. 

In summary, we have presented predictions for the in- 
clusive DPE heavy flavor production in heavy ion colli- 
sions at the LHC. The cross sections are large enough and 
can be investigated experimentally. The hard diffraction 
factorization was considered further corrected by absorp- 
tion corrections given by gap survival probability factor. 
For the Pomeron structure function, the HI diffractive 
parton density functions were considered. We investi- 
gate the theoretical uncertainty on the multiple interac- 
tion corrections for the nuclear case and addressed the 
diffractive coherent scattering. The main results are the 
estimations for PbPb collisions, where it has been ob- 
tained = 32.5 (0.32) /ib for charm (bottom). We 
verified that the diffractive channel dominates over ex- 
clusive photoproduction channel in proton-proton case, 
whereas it is two or more orders of magnitude smaller 
than the photon-Pomeron contribution in heavy-ion col- 
lisions. The enhancement of ultra-peripheral collisions 
compared to single diffractive channel is driven by their 
distinct A-dependences in each process. The two-photon 
contribution is also important and it is clear that an ac- 
curate study using relevant kinematic cuts (rapidity gap 
separation and transverse momentum spectrum) in the 
LHC is in order. 
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